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ABSTRACT 


1.  This  rsport  deals  with  Che  problems  of  variation  of  Pn  amplitudes 
in  the  regional  and  near- regional  distance  ranges  (200^2100  Kilo¬ 
meters),  The  data  used  were  recorded  by  Long  Range  Seismic  Measurement 
(LRSM)  vans  cf  the  VELA  Seismological  Center  as  a  result  of  earth¬ 
quakes  throughout  the  United  States  and  numerous  nuclear  and  chemical 
explosions  in  the  same  region. 

2.  It  is  shown  that  the  patterns  of  Pn  amplitudes  versus  L  in  Western 
and  Eastern  United  States  (UUS  and  EUS)  are  markedly  different  and  that 
these  differences  are  related  to  different  velocity  structures  in  the 
two  regions.  These  differences  extend  to  at  least  ISO  kilometers  depth. 

3.  Neither  the  WUS  nor  EUS  patterns  conform  even  approximately  to  that 
predicted  or  suggested  by  Gutenberg  and  Richter  in  the  1000-2000  kilo¬ 
meter  range.  The  hypothesized  shadow-zone  does  not  exist  and  over¬ 
estimation  of  magnitude  by  as  much  as  1.5  magnitude  units  is  frequently 
done  because  of  failure  to  properly  understand  the  patterns  of  radiation. 

4.  By  proper  calibration  of  WUS  by  use  of  numerous  events,  it  is  now 
possible  to  get  consistent  estimates  of  magnitude  at  all  distance  ranges 
for  most  explosions  and  earthquakes.  This  is  very  Important  when  body 
wave  magnitudes  are  used  as  an  essential  element  in  an  identification 
cri terion. 

5.  ODtoining  consistent  estimates  of  amplitude  as  a  function  of 
distance  for  a  particular  event  requires  a  knowledge  of  the  energy 
distribution  between  the  several  refracted  phases  ured  between  200  and 
2500  kilometers  distance.  Data  on  hand  show  that  the  energy  partition 
function  is  reasonably  uniform  throughout  the  regions  investigated  but 
that  locally  ir  may  vary  radically,  resulting  in  a  ten-fold  change  in 
relative  exclcatlon  of  two  refracted  phases. 

6.  Patterns  of  energy  radiation,  such  as  those  for  EUS  and  WUS,  are 
probably  approximately  determinable  from  knowledge  of  velocity  structure 
and  vice  versa  but  are  not  determinable  from  earthquake  Investigations. 
They  are  only  determinable  from  explosion  data  where  origin  time  and 
focus  are  accurately  known. 
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INTRODUCTION 


1.  Little  advance  has  occurred  In  our  understanding  of  so-called 
magnitude  values  computed  at  regional  or  near-regional  distances 
ov  use  oi  d y  waves  •'into  the  work  of  Richter  and  Gutenberg  (VESIAC 
Report  4alO-71-X,  Proceedings  of  the  VESIAC  Conference  on  Seismic 
Event  Magnitude  Determination,  May  1964).  Work  done  by  Romney  some 
time  ago  and  described  In  detail  In  this  report  constitutes  one  of 
the  few  exceptions  to  this  generalization.  Seismologists  have  tended 
to  simply  adopt  the  curves  given  by  Richter,  extend  them  by  arbitrary 
assumptions,  or  to  ignore  the  whole  matter.  The  most  cursory  glance 
at  available  Ampl ltude/Period  (A/T)  versus  Epicentral  Distance  (A) 
data  impresses  one  with  the  profound  control  exerted  on  such  data  by 
regional  crustal  characteristics.  Imposition  of  a  '‘standard"  A/T 
vi  rsus  A  c  hi  ve  on  dam  Tor  events  east  of  the  Rocky  Mountains,  here- 
alt*r  rcierrid  in  .is  Eastern  United  States  (EUS),  with  the  associated 
plotting  oi  magnitude  (m)  versus  A  rather  than  A/T  versus  A  obscured  the 
simple  A/T  versus  A  relationship  of  that  region.  Even  in  the  Western 
United  States  (WUS)  the  standardized  approach  to  magnitude  computa¬ 
tions  as  exemplified  by  VELA  reports  dealing  with  individual  Nevada 
Test  Site  (NTS)  events  (so-called  "shot  reports")  effectively 
obscured  significant  aspects  of  travel-time  and  amplitude  data 
as  well  as  evidence  for  the  vertical  velocity  structure  in  the  crust 
and  upper  mantle.  The  whole  concept  of  the  "crust"  of  the  earth  as 
defined  by  seismic  velocity  data  may  require  revision  based  on  data 
obtained  by  the  VELA  program.  This  report  attempts  to  advance  under¬ 
standing  o i  th<  problems  associated  with  magnitude  determinations  at 
regional  and  near • regional  distances.  Corollary  problems  that  became 
evident  during  the  analysis  will  be  discussed  at  the  appropriate 
places- 

2  Roth  short  and  l-'n;,  period  amplitude  measurements  of  high  quality 
for  both  earthquakes  and  explosions  arc-  now  available  as  a  result  of 
the  ARPA/AFIAl  program  They  show  that  the  A/T  versus  A  relationship  in 
the  range  200-2100  kilometers  tor  the  first  P  phase  as  published  in 
tabulated  toim  in  the  VELA  shot  reports  represents  a  comparatively 
crude  jt'empr  to  express  the  observed  pattern  of  variation  of  P 
amplitudes  with  distance.  That  curve  bears  no  relation  to  patterns 
of  P  energy  resulting  from  events  in  the  EUS.  The  supposed  profound 
decrease  and  subsequent-  increase  of  amplitude  in  the  neighborhood  of 
1000  kilometers  cpicentral  distance  (the  "observational"  basis  for  the 
"Shadow  Zone")  in  the  WUS  has  little  basis  in  fact.  The  misconception 
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arose  by  unknowingly  attempting  to  analyze  the  A/T  versus  A  data  for 
three  independently  propagating  refracted  phases  by  one  curve.  It 
will  in  shown  ui.»c  magnitudes  computed  according  10  the  formula 
appropriate  to  each  phase  largely  remove  the  apparently  profound 
dependence  of  determined  magnitude  upon  distance  of  observation. 

Systematic  variations  of  1.9  magnitude  units  occur  when  using  the 
shot  report  A/T  versus  A  table  at  regional  and  near- regional  distances 
(for  example,  see  Figures  8,  9,  10,  58,  60,  62,  and  all  plots  of  data 
from  earthquakes).  Such  overestimation  of  magnitude  becomes  a  matter 
of  concern  to  bomb  detection  and  identification  studies  as  some  pro¬ 
posed  identification  criteria  demand  a  valid  and  meaningful  ordering 
of  seismic  events  on  a  body  wave  magnitude  basis.  For  example,  mis¬ 
calculation  of  |he  magnitude  of  events  can  effectively  garble  such  a 
criterion  as  ARl '  versus  body  wave  magnitude.  Most  of  this  report 
deals  with  the  clarification  of  the  problem  of  understanding  the 
variation  of  P  amplitudes  as  a  function  of  distance  in  the  regional 
and  near-regional  distance  ranges  (150-2100  kilometers). 

3.  The  data  of  this  study  refute  widely  held  ideas  about  the  functional 
relationship  between  amplitudes  of  P  phases  and  epicentral  distance.  In 
so  doing,  the  data  clearly  show  that  some  current  ideas  about  the  proper 
placement  of  seismometers  for  optimum  detection  capability  are  grossly 
misleading  for  the  WUS  and  are  untrue  for  the  EUS  and,  by  extension,  for 
much  of  the  world.  The  WUS  represents  only  one  of  many  possible  areas  in 
terms  of  velocity  structure  and,  unfortunately,  most  studies  of  amplitude 
as  a  function  of  distance  have  treated  WUS  events.  Misinterpretation  of 
the  data  resulted  in  invalid  conclusions  for  even  that  region.  We  will 
demonstrate  the  complete  irrelevancy  of  the  currently  accepted  A/T  versus 
A  assumptions  to  conditions  in  EUS  though  many  seismologists  assume  these 
assumptions  to  be  valid  No  shadow  zone  exists  in  EUS  and  amplitudes  of 
the  Pa.  5  phase  increase  continuously  from  tcleseismic  distances  to  the 
epicenter  (that  is,  to  the  inner  limit  of  existence  of  the  phase  at  a 
hundred  or  so  kilometers).  A  10-fold  increase  in  signal-to-noise  ratio, 
relative  to  that  at  an  epicentral  distance  of  65°,  can  be  obtained  by 
observing  Pio.A  at  20°  or  Pg  5  at  11°. 
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U.  The*  observed  1 1  a>--c  1 -lime  c ompl i eat  ions  in  Wl’S  at  distance,  o:  It 
than  2000  t-  1  lometi  ■:  .  aru  largely  nonejvi  stent  :or  many  other  areas,  the 
general  belief  in  the  universality  of  such  corr.pl  1  cat  ion-,  results  t  rom 
confusing  the  pro'  lem  of  multiple  refractors  with  that  of  regional  varia¬ 
tions  of  refraction  velocities.  Possibly,  highly  seismic  regions  have, 
as  a  corollary  feature,  complex  velocity  structures. 

5.  This  repo-t  purports  to  demonstrate  the  general  statements  made  above 
ard  to  present  the  quantitative  relationships  between  amplitude,  distance, 
and  phase  for  FI'S  and  WUS.  We  will  consider  initially  El’S  because  of  its 
remarkable  simplicity  and  will  then  attempt  to  explain  the  major  problems 
inherent  in  understanding  amplitude  relationships  in  the  P  phases  arising 
from  WHS  events.  Any  careful  attempt  to  develop  a  meaningful  A/T  versus 
A  relationship  requires  the  ide~t it ication  of  phases,  the  determination 
of  thei'  distance  dependence-  and  their  average  energy  content  relative 
to  other  phases. 

6  Magnitude  has  varied  definitions  with  none  universally  accepted. 

Since  the  important  points  to  consider  include  relative  ordering  of 
events  and  establishment  of  their  relative  magnitude,  we  can  accept  any 
reasonably  meaningful  poi^t.  of  departure  for  our  definition.  The 
definition  of  the  magnitude  we  accept  is  as  it  is  deii-ed  by  several 
observations  at  20-30  ,  accordirg  to  Richter's^  published  curve  for 
P  body  waves.  Formulas  derived  herein  yield  a  uniform  estimate  of  the 
magnitude  by  use  of  data  obtained  between  200  kilometers  and  teleseismic 
distances.  We  can  then  use  those  formulas  to  determine  a  meaningful 
estimate  of  relative  magnitude  of  e  ents  over  any  distance  range. 

7.  I  he  suhs,  ripe  nonx  ->c  lat  jr  .■  used  to  characterize  particular  phases 
should  not  be  *.  al-vt  too  literally,  A  subscript  of  (7.9)  mea'vs  that  we 
attempt  to  treat  under  one  phase  and  equation  prof i It  s  with  observed 
velocities  m  the  ra^g«  C  /  —  7 . 9 > .  The  (8  5)  subscript  should  be  con¬ 
sidered  as  indication  of  velocities  in  the  range  (8.3-8  7)  while  (8.1) 
covers  (8.0-8. 1)  in  WUS.  The  P^o, r  s^bscr ipt  identities  the  first 
segment  oi  the  travel-time  curve  observed  beyond  Pg  5.  These  phases 
are  distinct  and  art-  associated  with  distinct  refractors. 

EASTF.RN  UNITED  STATES 

1.  Travel  -lime  Curves  for _ EUS.  The  standard  (T  -  A/8.1  versus  A)  fig¬ 

ures  of  the  shot  reports  show  clearly  the  essential  character  of  pn  travel¬ 
time  curves  for  EUS  (Figures  2  through  7).  The  Jef f reys-Bul len  (J-B) 
travel-  time  cic  'o  is  included  in  all  figures.  Figures  2  through  7  give 


•*See,  for  example.  Rivhtei,  C.  F.,  i 958,  Elementary  Seismology,  Freeman 
and  Company,  San  Francisco,  p.  688. 
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the  first  arrival  data  observed  for  the  indicated  events,  numerous  figures 
being  presented  in  order  to  make  it  clear  that  this  pattern  of  travel-time 
profile  typifies  epicenters  throughout  CUS.  Figure  1  indicates  the  epi¬ 
centers  of  these  events  and  all  others  considered  in  this  report.  Note 
that  the  events  ct  EUS  are  distributed  from  Eastern  New  Mexico  to  the 
Western  Atlantic  Ocean  and  from  the  South  Dakota-Nebraska  Border  to  the 
Louisiana  Gulf  Coast.  Clearly,  the  simple  profile  of  a  Pn  leg  of  approxi¬ 
mately  8.5  velocity  extending  from  around  150  to  about  2100  kilometers 
typifies  the  region.  At  that  distance,  a  deeper  refraction  becomes  first 
arrival,  travel  times  then  following  the  J-B  curve  with  regional 
characteristics  no  longer  playing  a  significant  role.  A  comparison  of 
this  ttav<‘1  time  curve  pattern  with  the  incomplete  pattern  of  curves  of 
WUS  presented  in  Figure  23  indicates  the  gross  contrast  of  velocity 
structures  of  the  two  areas.  More  discussion  of  this  point  will  follow 
a  discussion  of  the  WUS  profile.  For  the  moment,  we  note  only  the 
existence  of  tin  contrast  with  its  implications  of  different  patterns 
of  amplitudes  as  a  function  of  distance  from  epicenters  in  the  two  areas. 

2 .  .-l.-gm  tudes  as  computed  by  the  Shot  Report  Table  for  EUS  Events. 

The  aw  If  circles  on  Figures  8  through  16  indicate  magnitude  values  (msr) 
based  upon  the  standard  shot  report  A/T  versus  A  table  in  the  range 
of  1  2100  kilometers.  Clearly,  the  magnitude  estimates  are  a  function 

cf  distance  and  thu«  unsatisfactory  in  estimating  the  magnitude  of  EUS 
events.  Not  surprisingly,  a  generalized  A/T  versus  A  curve  attempting 
to  express  the  comp l c  xi  t  ics  of  WUS  velocity  structures  fails  completely 
vlu  n  implicit  t'  fUS,  a  simple  region  displaying,  essential  phase  con¬ 
tinuity  over  t  !u  omre  interval,  150-2100  kilometers. 

Ill  n<  4-  c  l  A/ r  on  o  and  the  Computation  of  Magnitude  for  El’S 

Fronts  . 

a.  It  is  a'-surud  that  the  relationship  between  m,  A/T,  and  A  is  of 
t  h*  term 


m  •  a  ♦  log  (A/T)  +  n  log  A  (1) 

where  m  is  magnitude,  a  and  n  are  constants  to  be  determined,  A/T  is 
measured  in  mi  I \  m  crons  per  second,  and  A,  the  epicentral  distance,  is 
measured  in  kilometer*  The  equation  can  be  rewritten  as 
t  -  m  •  a  •  leg  (A/T)  +  n  log  A  ,  the  correct  value  of  n  being  that  for 
which  (  is  independent  of  A  After  n  is  determined,  a  is  evaluated  by 

use  ot  the  nagii’tud*  value  determined  at  tcleselsmic  distances  (beyond 
th*-  Pg  5  ls>w  <f  tra.d  time  curve  or,  in  this  case,  beyond  2100  kilo- 
met  f  rs ) 

b  An  estimate  c.  t  n  can  be  obtained  by  determining  a  least  square  fit 
ot  log  A/T  vii'u.s  it  log  A  data  for  events  occurring  in  EUS.  The  six  events 
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selected  and  the  computed  values  of  n  with  their  standard  deviations  are: 


West  Virginia  Earthquakes  (25  Nov  1964) 

n 

-1.84  ♦  0.23 

Texas-Louisiana  Earthquake  (24  Apr  1964) 

-2.26  ♦  0.33 

GNOME  •  Eastern  Profile 

-2.07  ♦  0.39 

SALMON 

-2.37  ♦  0.34 

SS  Village 

-1.56  ♦  0.36 

Poplar  Bluff  Earthquake  (3  Mar  1963) 

-1.88  ♦  0.54 

Average  -2.0 

As  examples  of  the  influence  of  the  value  of  n  on  data  interpretation,  a 
detailed  presentation  of  the  data  for  the  first  two  events  above  is  given 
in  Table  I.  Note  that,  according  to  VELA  custom,  all  standard  deviations 
indicated  are  standard  deviations  of  a  single  determination.  Standard 
deviations  of  the  mean  are  generally  one-third  to  one-quarter  of  the 
standard  deviation  of  a  single  determination. 

TABLE  I 


West  Virginia  Earthquake  25  Nov  64  02:50:07.1  mts"  -  0-20 


Station  4  A/T 

(km)  nvt/s 

log 

4 

log  A/T 

log  A/T 
♦  log  4 

log  A/T  log  A/T  m 
+  2  logA  ♦  3  log  6 

m 

sr 

BRPA 

377 

51.8 

2.576 

1.714 

4.290 

6.866 

9.442 

3.60 

4.61 

CPSO 

403 

55.6 

2.605 

1.745 

4.350 

6.955 

9.560 

3.69 

4.75 

HDPA 

537 

14.5 

2.730 

.161 

3.891 

6.621 

9.351 

3.35 

4.64 

DHNY 

795 

19.7 

2.900 

1.294 

4.194 

7.094 

9.994 

3.82 

5.09 

LSNH 

1122 

3.26 

3.050 

0.513 

3.563 

6.613 

9.663 

3.34 

4.81 

GPMN 

1507 

2.93 

3.178 

0.467 

3.645 

6.823 

10.001 

3.55 

4.27 

WMSO 

1556 

4.36 

3.192 

0.639 

3.831 

7.023 

10.215 

3.75 

4.24 

RKON 

1784 

2.87 

3.251 

0.458 

3.719 

6.970 

10.221 

3.70 

3.36 

RTNM 

2022 

2.79 

3.306 

0.433 

3.739 

7.045 

10.351 

3.78 

3.35 

6.890 

3.62 

C  -  6.890  -  3.58  - 

a,  a  - 

-  3.31 
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Texas -louis  i ana  Earthquake  24  Apr  1964  07:33:53.5  mts  •  3.85  ♦  0.34 


Stat ion 

4 

(km) 

A/T 

log  6 

log  A/T 

log  A/T 
♦  log  4 

log  A/T 
+  2  log4 

log  A/T 
♦  3  log 4 

m8. 5 

msr 

JELA 

!  /() 

492.2 

2.2  30 

'2.692 

4,922 

7.152 

9.382 

3.88 

4.89 

GVTX 

333 

3/4.0 

2.522 

2.573 

5.095 

7.617 

10.139 

4.34 

5.27 

WMSO 

565 

14.1 

2.752 

1.149 

3.901 

6.653 

9.405 

3.38 

4.55 

EUAL 

574 

40.5 

2.759 

1.607 

4.366 

7.125 

9.884 

3.86 

5.01 

JUTX 

715 

52.0 

2,854 

1.716 

4.570 

7.424 

10.278 

4.15 

5.41 

CPSO 

883 

23.8 

2.946 

1.377 

4.323 

7.269 

9.466 

4.00 

5.38 

RTNM 

1127 

(2.04) 

3.052 

0.310 

3.362 

6.414 

9.980 

(3.14) 

4.61 

LCNM 

1212 

5.35 

3.084 

0.728 

3.812  • 

6.896 

10.249 

3.63 

4.93 

BLWV 

1  3  33 

7.49 

3.125 

0.874 

3,999 

7,124 

10,237 

3,85 

4.98 

BPPA 

16  34 

3-96 

3.213 

0.598 

3.811 

7.024 

7.070 

3.75 

3.80 

4.00 

C  -  7.070  -  3.85  -  a,  a  -  -3.22 

c.  Columns  6,  7,  and  8  give  values  of  C  for  n  ■  1,  2,  and  3, 

respectively.  The  values  of  Column  7  (n  ■  2)  appear  independent  of  A 

while  those  of  Columns  6  and  8  (n  ■  1  and  n  ■  3)  show  a  dependence  of  C 

upon  6  .  The  values  of  "a"  follow  immediately  (-3.31  and  -3.22) ,  their 
average  value  of  -3.27  being  adopted.  Finally,  the  formula  desired  reads 

meus  "  m8  5  "  *3’27  ♦  A/T  +  2  log  6  .  (2) 

This  equation  is  shown  graphically  in  Figure  100  as  the  Pg  5  curves,  A/T 
being  measured  in  millimicrons  per  second  and  A  in  kilometers.  Column  9 
gives  the  magnitude  computed  according  ro  this  formula  as  a  function  of 
distance.  Compare  these  values  with  those  of  Column  10  which  were  computed 

according  to  the  shot  report  table.  There  is  obvious  disagreement  of  msr 

values  with  the  teleselsmic  determination  (mrs)  while  the  mg  ^  values  yield 
the  statistically  Identical  estimate  of  magnitude  at  both  distance  ranges. 
Note  that  magnitude  estimates  that  differed  bv  1.74  units  (DHNY  and  RTNM) 
according  to  m&r  differ  by  only  0,04  units  when  computed  by  mg  5. 

d.  The  application  of  formula  (2)  to  other  events  in  EUS  achieves 
clarification  and  greater  internal  consistency  of  computed  magnitudes. 
Figures  8  through  16  give  the  relevant  data  for  all  events  studied.  The 
events  studied  include  surface  explosions  of  zero  focal  depth,  underwater 
explosions,  ard  earthquakes  of  unknown  but  shallow  focal  depth  The 
variation  of  typi  of  event  is  such  that  dependence  of  the  observed  patterns 
upon  a  specific  source  type  cannot  be  seriously  considered.  Also,  A/T 
appears  little  dependent  on  profile  azimuth.  The  only  GNOME  data  used  in 
this  part  of  the  report  comprise  the  profile  laid  out  eastward  from  the 
test  site.  Tne  only  event  showing  marked  discrepancy  between  mg, 5  and 

mtg  is  the  Lake  Superior  10-ton  shot,  mgj  being  4.3  ♦  0.35,  whereas  mts 
is  3.5  ♦  0.35.  The  reason  flor  this  discrepancy  cannot  be  demonstrated  but 
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it  will  be  shown  under  the  discussion  of  the  Nevada  Test  Site  (NTS) 
events  that  variation  in  rock  type  at  the  source  can  lead  to  large 
variations  in  the  pattern  of  energy  distribution  between  generated 
P-wave  signals. 

* 

e.  Note  the  systematic  discrepancy  of  the  msr  values,  maximum 
predicted  values  by  this  formula  being  around  1000  kilometers. 

Figure  100  shows  that  msr  overestimates  the  magnitude  by  1.5  units  at 
1000  kilometers  and  underestimates  it  by  0.4  units  at  1900  kilometers, 
a  difference  of  1.9  units  introduced  not  by  the  Jata  but  by  the  use  of 
an  inappropriate  formula. 

f.  Note  also  that  the  A/T  versus  A  data  of  these  eventfc  indicate 
the  total  lack  of  any  obvious  shadow  zone  phenomenon.  No  shadow  zone 
exists  in  EUS  and  the  supposed  problems  of  observing  at  distances  of 
less  than  2000  kilometers  because  of  low  signal  levels  do  not  exist  for 
EUS  events.  The  signal  level  at  1000  kilometers  attains  30  times  that 
predicted  by  the  msr  table  and  more  than  10  times  that  predicted  by 
msr  at  all  distances  between  400  and  1400  kilometers.  The  implications 
for  the  relative  advantage  of  regional  and  near-regional  stations  versus 
teleseismic  stations  for  regions  of  the  crustal  type  of  EUS  appear 
obvious. 
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g.  The  Vermont  Earthquake  data  depicted  in  Figure  16  provide  an 
excellent  example  of  the  misconceptions  that  can  be  reached  by  use  of 
inappropriate  formulas*  Very  few  stations  recorded  this  event , 
no  station  beyond  2017  kilometers  receiving  a  detectable  signal*  How¬ 
ever,  by  use  of  mgr  values,  one  obtains  an  estimated  magnitude  of  4.7, 
an  average  of  S  values  ranging  from  5.4  to  3.3.  An  estimate  of  4.7 
magnitude  on  the  basis  of  observations  at  regional  and  near-regional 
distances  is  in  pronounced  conflict  with  the  fact  of  total  lack  of 
recordings  at  teleseismlc  distances.  Such  an  event  should  yield  an 
A/T  value  of  10  nv  at  50°  and  should  be  widely  recorded  at  teleseismlc 
distances.  The  significance  of  this  misevaluation  will  be  pointed 
out  in  the  forthcoming  discussion  of  AR  versus  m.  Here  simply  note 
that,  by  use  of'ms.5  values,  the  spread  of  m  values  for  signals  clear 
in  both  period  and  amplitude  is  0.1  magnitude  units  as  compared  with 
2.1  for  msr  values.  The  estimated  magnitude  becomes  3.8  or  3.9, 
consistent  with  no  recordings  beyond  2000  kilometers. 

WESTERN  UNITED  STATES 

1.  Travel-Time  Curves  for  WUS.  In  the  WUS,  the  situation  becomes  much 
more  complicated  than  in  EUS.  Several  events  must  be  studied  in  order 
to  appreciate  the  complexity  even  to  the  extent  that  it  is  now  under¬ 
stood.  One  could  not  find  a  suite  of  events  which  would  establish  the 
basic  uniformity  of  the  WUS  as  such  uniformity  does  not  exist.  It  seems 
only  possible  to  appreciate  the  major  patterns  resulting  from  explo¬ 
sions  at  NTS  and  to  then  use  the  limited  data  available  from  other  areas 
to  suggest  the  similarities  and  dissimilarities  with  NTS  data  caused  by 
varying  eplcentrpl  positions  and  varying  azimuths  of  profiles.  We  shall 
consider  GNOME,  CLIMAX  (the  0.2  KT  chemical  explosion  at  Climax, 
Colorado),  50  shots  at  NTS,  and  17  earthquakes  distributed  from 
Central  Mexico  through  WUS  to  Queen  Charlotte  Islands. 


a.  It  is  convenient  to  begin  the  discussion  with  FORE  because  of 
the  excellent  and  significant  profiles  that  were  installed  for  this 
event.  See  Figure  17.  The  generally  observed  P79  refracted  wave  was 
well  obseived  from  200  kilometers  out  to  1000  kilometers  to  the  south¬ 
east,  reaching  Las  Cruces,  New  Mexico  (LCNM).  Remember  that  all 
arrivals  toted  on  figures  of  the  type  of  Figure  17  are  the  first 
observed  arrivals  at  the  several  stations.  No  observed  second  arrivals 
are  shown.  A  Pg  5  leg. was  received  as  first  arrival  at  Uinta  Basin 

Sei smological  Observatory  (UBS0),  Blue  Mountains  Seismological 
Observatory  (BMS0),  Forsyth,  Montana  (FRMA),  and  at  Glendive,  Montana 
(GIMA),  the  positioning  of  Raton.,  New  Mexico  (RTNM),  on  this  line  being 
probably  fortuitous.  Note  that  the  line  through  these  points  appears 
parallel  10  the  SS  Village  profile  but  with  a  time  intercept  on  the  zero 
distance  axis  of  14  seconds  rather  than  the  9  seconds  of  the  SS  Village 
profile.  The  time  arrivals  along  the  profile  to  the  ESE  provide  both 
complex  and  interesting  data.  Blanding,  Utah  (BXUT),  and  Durango, 
Colorado  (DRCO),  lie  well  below  the  two  pteviously  mentioned  refractions 
and  cannot  be  related  to  either  of  them.  More  discussion  of  these 
arrivals  will  be  presented  later.  RTNM  probably  did  not  actually 
receive  a  Pg  5  arrival  consistent  with  the  8.5  profile  through  FRMA  and 
GIMA  but  we  cannot  identity  the  phase  measured  at  this  time.  The  next 
stations  on  this  profile  appear  to  have  recorded  a  signal  with  a  slope 
of  about  8.5  but  with  a  time  intercept  on  the  zero  distance  axis  of 
nearly  20  seconds.  At  least,  they  did  not  record  the  phase  seen  at  BXUT 
and  DRCO  but  something  much  later  and  with  a  higher  refracted  velocity. 

b.  Finally,  the  P  phase  generally  received  as  first  arrival  beyond 
2000  kilometers  is  well  observed  at  EBMT  and  RKON,  2148  and  2338  kilo¬ 
meters  respectively,  as  well  as  at  Hempstead,  Texas  (HETX),  Hannah, 

North  Dakota  (HHND),  and  Ryder,  North  Dakota  (RYND),  the  latter  being  at 
an  epicentral  distance  of  1700  kilometers.  Thus  RYND  failed  to  observe 
earlier  but  less  energetic  ref factions,  recording  instead  the  strong 
PlO  5  arrival  a^  first  arrival. 

3.  Travel-Tim*  Curves  of  MISSISSIPPI.  MISSISSIPPI  (Figure  18)  shows 
similarity  to  FORE,  P -j  ^  being  observed  towards  LCNM  to  1200  kilometers 
and  to  shorter  distances  in  other  azimuths.  Hailey,  Idaho  (HLID),  and 
DRCO  received  as  first  arrivals  signals  preceding  the  expected  P7  9 
event.  The  profile  across  Colorado  and  the  Dakotas  follows  the 
previously  mentioned  8.5  leg  to  WNSD  but  then  jjmps  to  Pjq.5  ^or 
Academy,  South  Dakota  (AYSD),  at  1601  kilometers.  Other  points  are 
associated,  in  some  cases,  with  unclarified  phases. 
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4.  Travel -Time  Curves  of  BILBY  and  CLEARWATER.  The  interpretation  of  the 
BXUTt  HLID,  and  DRCO  arrivals  becomes  clear  as  soon  as  larger  events  such 
as  BILBV  and  CLEARWATER  are  considered  (Figures  19  and  20).  Thess  events 
show  the  existence  of  an  8.0  -  8.1  refraction  extending  eastward  from  NTS 
through  BXUT,  DRCO,  TDNM,  and  RTNM.  Shamrock,  Texas  (SKTX),  did  not  record 
this  phase,  a  more  energetic  phase  arriving  at  essentially  the  predicted 
time,  while  the  figures  show  clearly  that  Apache,  Oklahoma  (APOK), 
Grapevine,  Texas  (GVTX),  and  Durant,  Oklahoma  (DUOK),  observed  other  phases 
as  first  arrivals.  As  will  be  noted  below,  much  of  the  supposed  difference 
between  MISSISSIPPI  and  BILBY  and  the  inability  of  a  formula  based  upon 
MISSISSIPPI  data  to  reasonably  interpret  BILBY  data  results  from  different 
patterns  of  observation  in  a  heterogeneous  region. 

5.  Regional  Variations  of  Travel-Time  Curves  in  WUS.  Regional  variations 
within  the  WUS  appear  when  events  from  other  than  NTS  sites  are  considered. 
For  CLIMAX  (Figure  21),  located  in  Colorado,  the  first  arrivals  at  most 
stations  in  WUS  were  an  8.1  refraction  with  no  stations  measuring  a  7.9 
event  as  first  arrival.  As  noted  above,  the  eastern  profile  from  GNOME 
(Eastern  New  Mexico)  recorded  only  P8.5  while  Pg  ^  and  P79  were  seen  to 
the  west  and  northwest  (Figure  22),  the  Ps.l  event  being  observed  only 
very  near  GNOME.  Note  that  the  intersection  of  the  Ps.5  segment  through 
Montana  with  PlO.5  occurs  at  approximately  2000  kilometers  while  for  the 

SS  Village  profile  these  two  refractions  intersect  at  about  2200  kilo¬ 
meters.  The  suggested  P8.5  segment  through  Texas  and  Okalhoma  for  NTS 
events  intersects  Pl0.5  at  about  1700  kilometers.  P8.5  and  P7.9  inter¬ 
sect  at  a  distance  of  1250  kilometers  to  the  southeast  and  approximately 
1000  kilometers  to  the  north  and  northeast. 

6.  Conclusion.  Thus,  one  may  observe  at  least  six  different  P  phases 
from  NTS  shots  and  shots  from  different  localities  give  rise  to  different 
patterns  of  travel-time  curves.  Corppare  this  with  EUS  where  essentially 
one  travel-time  curve  for  one  phase  satisfied  all  observed  profiles.  Not 
surprisingly,  amplitude  interpretation  which  does  not  depend  upon  phase 
identification,  but  which  uses  a  single  A/T  versus  A  on  such  complex  data 
as  that  recorded  from  NTS  and  other  WUS  foci  has  yielded  nearly  chaotic 
results. 
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ESTABLISHMENT  OF  A/T  VERSUS  A  CURVES  FOR  WESTERN  UNITED  STATES 


I.  P7  o  Arrivals.  Romney  extensively  employed  P7 9  arrivals 
(VSC  Technical  Note  10)  to  determine  magnitudes  of  the  smaller  NTS 
events*  He  observed  this  phase  as  first  arrival  from  NTS  to  LCNM 
and  beyond  to  a  distance  of  approximately  1050  kilometers.  Northward 
from  NTS,  another  excellent  profile  extends  through  Mina,  Nevada 
(MNNV) ,  to  Pendleton,  Oregon  (PTOR),  including  HLID  for  events  in 
which  Pg.i  is  too  small  to  be  recorded.  PTOR  frequently  misses  P7(q 
in  favor  of  Pg  5.  Beyond  PTOR  In  Washington,  Montana,  and  the 
Dakotas,  either  Pgj  or  Piq.5  records  as  first  arrival  measured. 

DRCO,  to  the  east  of  NTS  never  records  P7,9  so  that  the  data  of  that 
station  should  be  ignored  when  considering  P7  9.  The  dependence  of 
A/T  on  A  for  the  NTS -LCNM  profile  (data  from  FISHER,  DORMOUSE,  STOAT, 
etc.),  the  LCNM-NTS  profile  (data  from  GNOME),  and  the  NTS-PTOR  pro¬ 
file  (numerous  events)  is  Identical  in  all  cases.  This  energy  path 
should  serve  as  an  excellent  source  of  internally  consistent 
estimates  of  magnitude  of  events  within  the  region.  The  procedure 
adopted  by  Romney  for  calibrating  this  energy  event  against  tele- 
seismic  data  is  as  follows  (Romney,  personal  communication): 

"a.  (U)  Since  1961  recordings  from  numerous  underground  nuclear 
explosions  at  the  NTS  have  been  made  at  mobile  stations  operated  under 
the  VELA  LRSM  program.  These  stations  have  been  located  at  many 
distances  and  azimuths  from  the  test  site.  Some  of  the  stations  have 
operated  more  or  less  continuously  and  have  recorded  data  from  many 
explosions,  while  others  have  moved  from  time  to  time,  recording  only 
a  limited  number  cf  explosions  at  a  given  location.  The  explosions 
have  included  yields  smaller  than  )  KT  and  larger  than  100  KT.  Shot 
points  have  been  somewhat  dispersed  within  the  test  site,  with  the 
result  that  the  propagation  paths  to  fixed  stations,  and  hence  the 
amplitudes,  have  been  affected  somewhat.  Variations  in  coupling 
which  appear  to  depend  upon  the  characteristics  of  the  rock  and 
other  conditions  near  the  explosions  have  also  been  noted.  It  Is, 
therefore,  necessary  to  correct  the  measurements  to  common  source 
conditions  (explosive  yield  and  coupling  factors)  before  the  data 
can  be  combined  to  give  a  good  representation  of  the  amplitude  versus 
distance  relationship. 

"b.  (U)  In  the  course  of  analyzing  recordings  from  the  explosions 
it  was  found  that,  for  any  given  station,  the  values  of  log  A/T  versus 
log  Y  (Y  -  yield)  fit  an  approximate  straight  line  if  the  measurements 
were  restricted  to  explosions  in  Yuc<a  Flats.  This  made  it  possible  to 
find  A/T  at  a  standard  yield  for  each  station  by  fitting  the  best 
straight  line  to  a  plot  of  amplitude  versus  yield.  The  amplitude  of 
10  KT  was  read  from  this  line  and  used  to  represent  the  best  value  of 
A/T  at  the  average  distance  of  the  station  from  t-'8  recorded  shots. 
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It  was  r^quirsd  ihat  at  least  lour  shots  which  bracket  10  KT  oe  recorded 
before  the  data  w^e  used,  This  procedure  tends  to  smooth  out  scatter  in 
th*  data  cau«°d  hv  small  variations  in  coupling  or  in  the  path  of  pro¬ 
pagation  and  by  noise. 

"  (U)  Ifn  average  values  tog  */T  at  10  KT  for  each  station  were 

plotted  as  a  function  of  distance  (figure  101).  These  points  detine  a 
line  given  by 


log  A/T  -  9,50  ♦  0  72  -  (3.04  ♦  0,26)  log  A  (3) 

where  A  is  measured  in  kilomet  T  in  seconds ,  and  A  in  millimicrons. 
Firs',  arrivals  at  stations  between  200  and  1050  kilometers  were  used. 

*  The  slope  of  this  curve  is  almost  identical  with  that  found  in  LOGAN 
and  BLANCA 

"d.  (IJj  The  amplitudes  ot  P_  _  can  he  associated  with  telescismic 
magnitudes  by  comparison  with  the  signals  recorded  at  long  ranges  from 
the  LRSM  stations.  The  basic  da  a  used  came  from  the  13  explosions 
which  were  recorded  at  five  or  mere  stations  at  teleseismic  distance 

.  i'j 

(  A>  16  )  as  well  as  in  the  P*i  range.  These  shots  v<^re  generally 
larger  than  10  KT  since  only  these  larger  shots  can  be  satisfactorily 
recorded  at  t*-leseismic  ranges  .  Unified  magnitudes  w°r..-  computed  for 
each  such  explosion  from  the  mcasur-.d  amplitudes  and  p-riod  of  P.  The 
average  magnitude  for  each  e\?rt  was  then  compared  with  the  value  of 
log  A/T  at  500  kilometers,  a?  found  by  a  least  square  fit.  of  the  P_  Q 
amplitude  versus  di>«ance  data  Results  are  shown  in  Figure  102  it 
may  be  seen  'hat  the  otagnirude  ;ncitase«  more  rapidly  than  log  A/T  at 
500  kilometer1,  The  reiatjon«hip  was  tound  by  least  squares  *o  be 

m?  9  -  2,36  ♦  0.2  5  t  (1.21  ♦  0 . 1 1 )  log  (A/T>500  (4) 

This  implies  that  the  '-■'le5°ism?c  amplitudes  increase  more  rapidly 
with  yield  than  does  the  amplitude  at  500  k i  Icon  t »*rs  This  probably 
results  from  the  shift  toward  l r require tes  n  the  source  as  the 
yield  incteas**  together  with  the  reduced  attenuation  at  low  fre¬ 
quencies. 

**«-.  (u)  Equations  (3)  and  (4)  mav  be  combined  to  give  the  P_ 
amplitudes  as  a  function  of  magnitude  The  result  is  7’ 

m?  9  -  -7,55  «■  1.21  (log  A/T  ♦  3,04  log  A  ),"  (5) 
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Note  on  Figure  100  that  a  discrepancy  between  msr  and  my  g  of  0.6  occurs 
at  my  g  equal  to  5.0,  my  g  being  smaller  than  msr. 

2.  £ 10»5  Arrivals.  We  shall  fit  a  smoothed  curve  through  the  Gutenberg - 

Richter  curve  between  2100  and  3000  kllpmeters  and,  in  addition,  obtain 
the  values  of  amplitude  for  Pio.5  observed  at  1600  kilometers  in 
North  Dakota  at  RYND,  etc.  The  equation  obtained  is 

•  -10.35  +  log  A/T  ♦  4  log  A  (6) 

where  A/T  is  measured  in  millimicrons  per  second  and  A  in  kilometers. 

The  corresponding  curve  is  shown  on  Figure  100  as  Pio.5*  As  far  as  is 
known,  this  dependence  of  A/T  upon  A  occurs  independently  of  azimuth. 

3*  £8.5  Arri  vals .  As  noted  above,  PTOR,  at  approximately  1000  kilo¬ 

meters,  generally  records  Pg5  os  fust  measured  arrival.  A  phase  with 
this  velocity  is  observed  as  first  arrival  in  Montana  and  the  Dakotas. 

To  the  southeast,  stations  in  Texas,  Oklahoma,  and  beyond  observe  Pg^ 
as  first  arrival  to  1800-2000  kilometers.  The  time  intercepts  for  these 
three  Pg  5  profiles  appear  quite  different,  both  from  each  other  and 
from  that  typical  of  this  refraction  in  EUS.  Even  so,  we  proceed  by 
attempting  to  employ  the  meus  formula  for  these  arrivals.  The  results 
agree  satisfactorily  with  the  magnitude  computed  for  P79  and  tele- 
seismic  distance.  Therefore,  we  adopt  the  meu5  (i.  e.,  m^  5)  curve  for 
computing  magnitudes  for  arrivals  on  Pg.5  profiles  in  WUS. 

4.  £8.1  Arrivals .  When  establishing  a  curve  for  mg  ^ ,  a  most  interest¬ 

ing  result  appears  in  this  case.  A  well-observed  profile  for  this 
refraction  for  NTS  shots  require*  a  large  explosion,  only  DRCO  recording 
this  phase  for  most  events.  Either  this  phase  possesses  very  little 
energy  or  the  energy  dissipates  rapidly  along  the  refractor.  In  this 
connection,  remember  that  for  the  CLIMAX  shot  in  Colorado,  a  small 
event  generated  a  widely  obser zed  Pg.j  refraction.  A/T  values  from 
CLIMAX  measured  toward  Arizona  d-crcased  as  about  the  inverse  third 
power  of  the  dis’ance  while  along  the  profile  from  NTS  through 
Southern  Colorado  into  Northern  New  Mexico,  A/T  values  for  Pg.j  arrivals 
decrease  as  the  inverse  seventh  power  of  the  distance.  Note  that  DRCO 
occurs  op  both  the  Cl.TMAX-Arizona  and  the  NTS  New  Mexico  profiles.  High 
amplitudes  typifv  Pg  1  arrivals  at  short  distances  from  NTS.  The  signal 
decays  so  rapidly  with  Increasing  distance  that  it  disappears  into  the 
noise  beyond  DRCO  for  small  events.  When  one  attempts  to  derive 
equations  to  use  along  the  NTS-Ncw  Mexico  profile  by  use  of  CLEARWATER 
and  BII.BY,  it  becomes  apparent  that  the  data  of  the  two  events  require 
a  0.9  difference  in  the  leading  coefficient  of  the  derived  equations. 
Thus,  for  CLEARWATER,  shot  in  mesa  tuff,  or.e  obtains 

mg  1  (CLEARWATER)  -  - 1 7 . 20  +  log  A/T  +  7  log  A  (7) 
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while  for  BILBY,  shot  at  2314  foot  depth  in  valley  tuff, 

m_  .  (BILBY)  -  -16.30  ♦  log  A/T  ♦  7  log  6  .  (8) 

0  •  1 

Apparently,  CLEARWATER  coupled  eight  tinea  as  much  energy  (relative  to 
that  coupled  into  the  other  measured  refracted  phases)  into  this  phase 
as  did  BILBY.  If  one  computes  DRCO  P„  .  arrivals  for  events  with  vary¬ 
ing  shot  locations  according  to  Pg  ,  ^CLEARWATER)  and  Pg  ;  (BILBY),  the 
following  pattern  emerges: 


TABLE  IV 


Depth  . 

Magnitude 

DRCO  by  mg  ^ 
(CLEARWATER  5 

DRCO  by  i 
(BILBY)  . 

DORMOUSE 

Alluvium  (3) 

1191  ft 

3.64 

3.72 

MISSISSIPPI 

Valley  Tuff  (9) 

1615 

4.76 

- 

4.89 

CIMARRON 

Alluvium  (9) 

1000 

3.98 

3.71 

- 

STILLWATER 

Alluvium  (9) 

625 

3.46 

3.27 

- 

CODSAW 

Alluvium  (9) 

696 

3.40 

3.03 

m 

MINK 

Alluvium  (3) 

630 

2.72 

3.06 

- 

HAYMAKER 

Alluvium  (12) 

1351 

4.52 

- 

4.46 

DES  MOINES 

Mesa  Tuff  (12) 

3.72 

3.71 

m 

MARSHMALLOW 

Mesa  Tuff  (12) 

3.99 

3.63 

m 

MADISON 

Mesa  Tuff  (12) 

3.71 

3.73 

m 

Thus,  all  of  the  mesa  tuff  shots  behaved  similarly  to  CLEARWATER  but  so 
did  all  of  the  shallow  valley  shots.  Only  HAYMAKER  and  MISSISSIPPI 
coupled  as  little  energy  into  this  refractor  as  did  BILBY.  Thus,  the 
use  of  data  obtained  from  this  refraction  throughout  WUS  requires 
extensive  crustal  calibration  (CLIMAX  profile  versus  NTS  profile)  as 
well  as  calibration  of  energy  in  the  phase  relative  to  source  location. 

In  other  words,  this  phase  does  not  appear  very  useful  for  routine 
Investigations.  The  very  high  rate  of  amplitude  decrease  eastward  from 
NTS  together  with  the  markedly  different  behavior  on  another  azimuth 
through  the  same  area  imply  the  disappearance  eastward  of  the  velocity 
contrast  requisite  for  the  refraction.  Such  an  azimuth  dependent  phenomnon 
requires  lateral  crustal  inhomogeneities. 

DETAILED  ANALYSES  OF  STATION  MAGNITUDE  VALUES  FOR  MISSISSIPPI,  BILBY, 

AND  HARDHAT 

1.  Introduction.  Three  events  will  be  discussed  in  detail.  Two  of 
them,  MISSISSIPPI  and  BILBY,  illustrate  the  problem  of  observing  different 
profiles  and  thus  possibly  different  phases  for  different  events  while 
interpreting  all  data  by  an  inadequate  generalized  formula.  All  three 
events  illustrate  a  type  of  problem  that  even  proper  phase  assignment 
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will  noc  defeat.  Obtaining  the  same  estimated  magnitude  from  several 
phases  of  the  same  event  requires  that  the  energy  partition  function 
for  the  event  be  identical  to  that  for  the  event  or  events  used  in 

establishing  the  function.  Otherwise,  magnitude  becomes  a  function  of 

phase. 

r 

2.  MISSISSIPPI .  MISSISSIPPI  (Figures  18  and  64)  probably  represents  the 
most  significant  event  used  in  the  establishment  of  the  m^  „  formula. 

Not  surprisingly  then,  m^  ^  and  m  values  agree  closely,  use  of  the 

standard  EUS  m  ,  formula  obtaineaexcellent  agreement  with  m^.  The 
station-by-station  data  are  as  follows;: 

TABLE  V 

MISSISSIPPI 


A 

Magnitude 

m 

sr 

A 

Magnitude 

m 

sr 

MNNV 

234 

5.06 

(7.9) 

6.0 

SSTX 

1496 

4.33 

(8.5) 

5.0 

KNUT 

278 

4.68 

(7.9) 

5.5 

WNSD 

1507 

4.83 

(8.5) 

5.6 

FMUT 

408 

- 

? 

4.8 

HB0K 

1555 

4.98 

(8.5) 

5.5 

TFCL 

415 

3.70 

(7.9) 

4.5 

CKBC 

1578 

4.52 

(8.5) 

5.0 

WINY 

483 

4.83 

(7.9) 

5.8  * 

WMS0 

1596 

4.55 

(8.5) 

4.9 

FSAZ 

485 

4.86 

(7.9) 

5.5 

AYSD 

1601 

4.95 

(10.5) 

6.0 

CPCL 

490 

5.10 

(7.9) 

5.7 

HMBC 

1670 

4.84 

(8.5) 

5.0 

MVCL 

517 

- 

? 

•» 

PU0K 

1714 

? 

TFSO 

536 

4,67 

(7.9) 

5.0 

CT0K 

1910 

4.30 

(10.5?) 

4.3 

DRCO 

733 

4.89 

(8.1) 

5.5 

SJTX 

1970 

4.86 

(10.5) 

4.9 

HLID 

739 

4.61 

(8.1) 

4.6 

SEMN 

1971 

4.42 

(10.5) 

4.5 

BMSO 

863 

4.47 

(7.9) 

5,2 

MPAR 

2083 

4.49 

(10.5) 

4.6 

PTOR 

971 

4.88 

(8.5) 

5.3 

LR  17 

2099 

4.85 

(10.5) 

4.9 

LCNM 

1101 

4.41 

(7.9) 

5.2 

HTNM 

2120 

5.12 

(10.5) 

5.2 

PMWY 

1027 

4.40 

(7.9) 

5.2 

LR  16 

2460 

5.11 

(10.5) 

5.1 

HKWY 

1126 

4.51 

(8.5) 

6.0 

ARMS 

2460 

5.05 

(10.5) 

5.1 

MUWA 

1392 

4.27 

(8.5) 

5.3 

CVTN 

2565 

4.95 

(10.5) 

5.0 

Average  of  all  stations  at  less  than  1900  kilometers  -  4.76  (4. 3-5. 2) 
Average  of  28  stations  between  1900  and  4064  kilometers'  4.69  (4. 3-5.1) 
Average  of  all  m  values  at  less  than  1900  kilometers  -  5.3  (4. 3-6.0) 

The  figures  in  parentheses  are  the  extreme  values  in  each  group  of 
values. 
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Apparently,  computations  based  upon  phase  assignment  and  appropriate 
formulas  have  made  the  computed  magnitude  independent  of  epicentral 
distance  and  have  also  reduced  the  range  of  variation  of  the  data. 

Whereas  the  range  in  computed  values  was  4.3  to  6.0  using  the  shot 
report  table,  the  use  of  phase  curves  reduced  the  range  to  4.3  to  3.2, 
i.  e. ,  half  the  previous  range.  The  pairs  PMWY  -  HKWY  and  WMSO  -  AYSD 
clearly  show  the  influence  of  proper  phase  Identification.  For  the 
former  pair,  m  values  were  5.2  and  6.0,  while  properly  interpreting 
the  PMWY  arrival  as  P.  and  the  HKWY  arrival  as  Pg  ,  leads  to  magnitude 
estimates  of  m_  ■  4.40  and  mR  _  ■  4.51,  respectively.  The  m  values 
for  WMSO  and  AYSD  are  4.9  and  old,  while,  when  making  proper  p^ase 
identifications,  mfl  .  for  WMSO  Is  4.55  and  m.Q  .  for  AYSD  is  4.95. 
Regarding  the  rejection  of  such  stations  as  TFCl  and  MVCL,  we  argue  that, 
for  Interpretation,  we  require  a  reasonably  clear-cut  relationship 
between  data  obtained  at  a  given  station  and  that  recorded  at  tele- 
seismic  distances.  If  we  have  not  established  as  yet  such  a  relation¬ 
ship,  the  proper  mode  of  interpretation  of  the  data  of  such  stations 
is  unknown  and  we  ignore  their  recordings.  Further  investigations 
than  reported  here  may  render  the  data  of  such  stations  meaningful 
in  a  magnitude  sense,  and  thus  useable. 

3.  BILBY.  For  BILBY  (Figures  19  and  65),  the  data  are  as  follows: 

TABLE  VI 
BILBY 


A 

Magni tude 

m 

sr 

A 

Magni tude 

m 

sr 

MNNV 

242 

5.46  (7.9) 

6.0 

TDNM 

880 

5.59 

(8.1) 

5.3 

KNUT 

284 

5.42  (7.9) 

5.9 

RTNM 

1039 

5.32 

(8.1) 

5.2 

CPCL 

483 

?  early 

AZTX 

1278 

5.05 

(8.5) 

6.2 

WINV 

493 

5.01  (7.9) 

5.6 

FRMA 

1232 

5.17 

(8.5) 

6.3 

5.13  (8.1) 

TKWA 

1337 

7 

6.0 

MVCL 

520 

?  early 

SKTX 

1426 

5.26 

(8.5) 

6.2 

TFSO 

529 

5.68  (7.9) 

6.2 

GIMA 

1487 

5.05 

(8.5) 

5.9 

BXUT 

586 

5.68  (8.1) 

5.9 

WMSO 

1593 

5.18 

(8.5) 

5.5 

UBSO 

668 

5.46  (8.5) 

6.7 

RYND 

1705 

5.40 

(10.5) 

6.0 

DRCO 

732 

5.91  (8.1) 

5.9 

GVTX 

1794 

5.75 

(8.5)  1 

5.4 

HLID 

747 

5.51  (8.1) 

5.6 

DUOK 

1831 

5.81 

(8.5)  ? 

5.5 

BMSO 

871 

? 

4.24  (8.5)  ? 

5.6 

HHND 

1926 

5.83 

(10.5) 

6.0 

Average 

of 

18 

stations  at 

less  than 

1900  kilometers 

5.44  (5 

.1-5.9) 

Average 

i  of 

18 

stations  beyond  1900  kilometers 

(1926-9100)  - 

5.6  (5 

.2-6.0) 
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N  ir  e  that,  •  r  tr>-  '  :« ;  al  550  1040  kilometeis,  the  only  profile  observed 
for  BILhv  wa?  th-  MS  -  BXUT  -  DRCO  -  TDNM  -  RTNM  path,  and  that  Pg4l  was 
observed  as  tir>t  arrival  at  all  stations.  When  calculated  according  to 
n>7  9  (ra  itlv  ac?jming  that  all  arrivals  in  this  distance  range  are  £7.9)* 
magnitude  ««tlmat~s  at  these  stations  are  5.33,  4.64,  4,64,  and  4.55, 
respectively,  values  below  estimates  made  upon  P79  arrivals  and  values 
apparently  decreasing  with  distance.  The  values  ente-ed  In  the  tabulation 
above  w<  found,  as  described  previously,  by  calibration  against  the 
tele«eiemic  data  of  this  event.  Note  that  the  and  mg  5  values  appear 

svs' .mat ically  low  relative  to  m^Q  5  and  mts  values.  We  obtained,  by  use 
01  fh'.  same  formulas,  excellent  agreement  between  magnitude  estimates  at 
all  ranges  for  MISSISSIPPI.  This  Implies  that  the  energy  pertition  for 
these  two  events,  BILBY  and  MISSISSIPPI,  differed  markedly.  The  apparent 
agreement  of  mg  ^  and  mts  occurs  because  we  calibrated  the  mg  j  formula 
by  use  of  this  event.  Probably,  calibration  against  107  9  and  mg  5  would 
have  had  more  meaning,  making  all  of  the  shallow  refractions  less 
energetic  relative  to  mAQ  5  and  deeper  refractions  for  BILBY  than  for 
MISSISSIPPI.  Such  variation  In  the  energy  partition  function  requires 
either  nonradial  symmetry  in  the  source  radiation  function,  varying 
patterns  of  vertical  velocity  structure  in  the  two  shot  areas  sufficient 
to  give  the  observed  differences  of  energy  level  in  the  lunes  of  the 
focal  sphere  giving  rise  to  the  several  refractions,  or  frequency  and 
depth  dependent  attenuation  factors. 

4.  HARDHAT.  For  HARDHAT  (Figure  75),  the  data  are  as  follows! 

TABLE  VII 
HARDHAT 


MNNV 

228 

4.96 

(7.9) 

5.96 

ATNV 

2r) 

5.30 

(7.9) 

6,11 

KNur 

288 

?o  1  2 

(7.9) 

5,89 

INCL 

336 

4.09 

(7.9) 

4.92 

WMAZ 

399 

4  90 

(7.0) 

5.49 

FMUT 

403 

4,06 

(7,9) 

4.83 

WJNV 

474 

5.19 

(7.9) 

5.78 

FSAZ 

490 

5.18 

(7.9) 

5.80 

CPCL 

500 

?  early 

MVCL 

51? 

?  early 

SFAZ 

5h9 

5.63 

(7  9) 

6.13 

VNur 

669 

5. 19 

(7.9) 

VTOR 

•  6  88 

4.90 

(7.9) 

SVAZ 

711 

4,47 

(7.9) 

HUD 

729 

3  48 
4.08 

(8.1CL) 

(8. IB) 

4,34 

DRCO 

734 

4,32 

4.92 

(8.1CL) 

(8. IB') 

5.20 

MLNM 

781 

3.56 

(8.1CL) 

4.29 

TCNM 

902 

4.47 

(8.1CL) 

4,66 

5.34 

(8, IB’) 

4.86 

PTOR 

961 

4.44 

(8.5) 

4,86 

LCNM 

1017 

3.79 

(7,9) 

4.64 

RTNM 

1043 

Very 

late 

TRWA 

1070 

4.27 

(8,5) 

5.74 

EPTX 

1095 

3.46 

(8.5) 

4.94 

EFTX 

1209 

2.89 

(8.5*) 

4,21 

GNNM 

1244 

3,30 

(8.5*) 

4.56 

BMTX 

1324 

4,05 

(8.5') 

5.18 

SSTX 

1501 

3.60 

(8.5') 

4.32 

WNSD 

1504 

4.24 

(8.5') 

4.96 

HBOK 

1557 

4.18 

(8.5) 
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TABLE  VII  (Contd) 


A 

Magnitude 

msr 

A 

Magnitude 

msr 

WMSO 

1598 

? 

• 

MPAR 

2084 

3.90  (10.5) 

5.00 

LPTX 

1765 

? 

NGWS 

2500 

4.13  (10.5) 

4.13 

SEMN 

SJTX 

1967 

1975 

4.41  (8.5) 
4.09  (10.5) 

4.02 

5.37 

MMTN 

2728 

4.34  (10.5) 

4.37 

The  three  symbols  for  Pg#5 

arrivals 

identify  probable  8.5  refract 

ions  from 

different  profiles  with  different  intercept  times.  The  identification  of 
the  primed  events  as  P8.5  Is  tentative.  Uncertainty  exists  as  to  the 
appropriateness  of  using  mg#^  (CLEARWATER)  or  mg#j  (BILBY)  (prime  imply¬ 
ing  0.3  adjustment  of  BILBY  formula  based  upon  discussion  of  BILBY).  It 
the  computed  magnitudes  are  grouped  according  to  phase,  the  following 
tabulation  is  obtained: 


*7.9 

*8.1 

*8.5 

*8.5* 

3 

00 

• 

*10.5 

4.96 

3.48 

4.08 

4.44 

4.05 

3.46 

4.09 

5.30 

4.32 

4.92 

4.27 

3.60 

2.89 

3.90 

5.12 

3.56 

4.16 

4.24 

4.18 

3.30 

4.13 

4.09 

4.74 

5.35 

4.41 

4.34 

4.90 

4.06 

5.19 

5.18 
5.63 

5.19 
4.90 
4.47 
3.79 

The  two  arrivals  that  were  1  second  early  relative  to  P7  9  and  had  very 
low  computed  m7^9  magnitudes  (4.09  and  4.06),  together  with  the  unusually 
low  estimate  at  LCNM  (3.79)  introduce  an  abnormal  spread  into  the  m 7  9 
data.  If  one  ignores  these  arrivals,  the  average  m7#9  -  5.17,  if  they 
are  included,  the  average  m7>9  -  4.87.  The  mean  mg#;  (B')  -  4.62,  t be 
mean  mg^  •  4.34,  the  mean  m^g^  -  4.12.  with  the  mgj#  and  1113  51  values 
even  lower.  Obviously,  the  fo’-mulas  used  for  MISSISSIPPI  do  not  satisfy 
HARDILAT.  No  obvious  distance  dependence  appears  in  any  phase.  The  path 
length  differences  between  P7>9,  Pg#c,  and  Pi0.5  represent  small  frac¬ 
tions  of  the  total  path  lengths  so  that  one  cannot  explain  a  dependence 
of  magnitude  upon  T>hase  as  being  due  to  a  higher  frequency  signal  for 
HARDHAT  than  for  MISSISSIPPI  with  a  more  rapid  decay  of  signal  with 
distance.  The  lack  of  dependence  of  magnitude  upon  distance  in  each 
phase  also  suggests  this  conclusion.  Another  explanation  seems  required. 
The  decrease  in  ^10. 5  values  shown  above  from  those  in  the  shot  report 
largely  results  from  miscalculation  in  the  shot  report. 
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COMPARISON  OF  MAGNITUDES  COMPUTED  BY  msr  AND  THE  SEVERAL  CURVES  OF  THIS 
REPORT 

A  comparison  of  Che  magnitudes  computed  by  msr  and  the  several  curves 
.  established  here  clarifies  Che  significance  of  employing  Che  proper 
curve  (see  Figure  100).  AC  1900-2000  kilometers,  mar  underestimates 
relative  Co  ms. 5  by  0.4.  AC  1600  kilometers,  msr  underestimates  Che 
magnitude  by  0.3  relative  to  ms. 5  and  1.0  relative  to  mio.5.  Even  if 
using  the  new  curves,  an  error  of  0.7  will  be  made  if  one  mis identifies 
the  phase  measured.  At  1100  kilometers,  msr  overestimates  with  respect 
to  both  relevant  curves,  0.5  relative  to  017,9  and  1.5  relative  to  mg. 5. 
Again,  if  017,9  and  mg. 5  are  used  but  misidentif lcation  of  a  phase  occurs, 
an  error  of  1.0  will  be  made  in  magnitude,  PT0R  providing  a  perfect 
exanple  of  this  phenomenon.  This  station  routinely  sees  Ps.5  as  first 
arrival  for  NTS  events  though  being  only  970  kilometers  from  NTS.  Note 
the  gross  discrepancy  between  msr  and  ms. 5  at  all  distances  less  than 
1400  kilometers.  Fortuitously,  these  curves  cross  around  1700  kilometers 
due  to  the  attempt  to  get  by  some  means  from  something  approximating  a 
mean  between  m7.9  and  mg, 5  at  1000  to  mio.5  at  around  1800  kilometers. 
Between  1000  and  2000  kilometers  the  msr  and  117.9  curves  nearly  parallel 
one  another,  the  difference  between  the  two  curves  being  0.63  at 
msr  "  5  and  0.95  at  msr  ■  3.  We  included  mg.i  (CLEARWATER)  on  Figure  100 
to  demonstrate  the  grossly  different  behavior  of  this  phase  with  distance 
from  that  of  all  other  curves.  Inside  600  kilometers,  Pg.i  (CLEARWATER) 
appears  larger  than  either  P7.9  or  Pg.5  for  m  •  4,  while  it  decays  to 
1.5  ny.  at  1000  kilometers  for  m  -  4. 

RECALCULATION  OF  ALL  SHOTS  PREVIOUSLY  ANALYZED  IN  SHOT  REPORTS 

Subject  to  the  limitations  and  uncertainties  of  the  type  just  described, 
the  A/T  versus  6  versus  magnitude  data  of  all  previously  analyzed  events 
were  recalculated  based  upon  the  formulas  given  above.  As  none  of  them, 
other  than  SHOAL,  were  in  granite,  one  does  not  expect  problem  of  the 
type  met  in  analyzing  the  data  of  HARDHAT.  Figures  24  through  79  present 
the  recalculated  magnitude  values,  the  figures  for  selected  events 
(Figures  8,  9,  10,  58,  60,  62,  64,  65,  73,  76,  and  78)  including  shot 
report  values  of  magnitude  for  comparison  and  study  purposes.  The  mag¬ 
nitude  of  these  events  was  obtained  by  averaging  all  recomputed  values, 
the  symbol  n^'  being  adopted  for  this  average  value.  The  values  of  mb' 
are  those  used  throughout  the  following  discussion  and  tables.  It  is 
clear  from  the  indicated  figures  that  the  gross  dependence  of  magnitude 
upon  distance  which  was  present  in  shot  report  values  of  magnitude  has 
been  largely  eliminated  by  use  of  the  phase  formulas.  Small  events  only 
recorded  at  less  than  1500  kilometers  are  reduced  in  magnitude  (m^') 
from  msr  values  by  as  much  as  a  full  magnitude.  The  following  tabulation 
(TABLE  VIII)  gives  all  mb'  values  computed  by  procedures  of  this  report 
for  explosions,  together  with  mgr  values  and  values  derived  by  Romney 
by  use  of  A/T  values  at  500  kilometers  for  the  smaller  events. 
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TABLE  VIII 


Magnitude  Data  on  Selected  NTS  Shota 


SHOT  mu  msr 

Alluvium 


SHREW 

1.2 

2.3 

MAD 

2.37 

3.7 

MINK 

2.72 

2.93 

3.7 

CHINCHILLA  II 

2.73 

2.84 

3.8 

ROANOKE 

2.94 

3.21 

3.9 

CHINCHILLA 

3.16 

3.33 

4.2 

ALLEGHENY 

3.04 

3.31 

4.0 

KAWEAH 

3.15 

3.36 

4. 1 

STILLWATER 

3.46 

3.40 

• 

CODSAW 

3.40 

3.54 

4.3 

BOBAC 

3.39 

3.60 

4.2 

RINGTAIL 

3.21 

3.36 

4.1 

SACRAMENTO 

3.32 

3.40 

4.1 

STOAT 

3.38 

3.60 

4.5 

WICHITA 

3.47 

3.58 

4.5 

AGOUTI 

3.55 

3.83 

4.5 

SANTEE 

3.30 

3.56 

4.2 

ARMADILLO 

3.51 

3.78 

. 

PACKRAT 

3.69 

3.73 

4.2 

CASSELMAN 

3.64 

3.96 

4.6 

HYRAX 

3.43 

3.60 

4.4 

ACUSHI 

3.51 

3.75 

4.4 

DORMOUSE  PRIME 

3.73 

3.91 

4.6 

DORMOUSE 

3.64 

3.83 

4.5 

CIMARRON 

3.98 

4.38 

YORK 

3.69 

3.92 

4.4 

PEBA 

3.54 

3.78 

4.4 

PASSAIC 

3.69 

4.06 

4.3 

FISHER 

3.45 

3.68 

4.3 

MERRIMAC 

3.80 

3.87 

4.3 

HAYMAKER 

4.52 

4.5 

4.9 

SEDAN  (crater) 

4.12 

- 

4.7 
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TABLE  VIII  (Contd) 


SHOT 

Valley  Tuff 

mb» 

"to 

mar 

AARDVARK 

4.55 

4.7 

4.9 

STONES 

4.50 

4.6 

4.8 

FORE 

4.85 

4.9 

5.2 

MISSISSIPPI 

4.76 

4.7 

5.1 

BILBY 

5.61 

5.6 

5.8 

Mesa  Tuff 

FEATHER 

2.13 

_ 

3.2 

CHENA 

2.76 

m 

3.9 

PLATTE 

3.44 

- 

- 

ANTLER 

4.03 

- 

4.7 

MADISON 

3.71 

- 

4.5 

DES  MOINES 

3.72 

- 

4.5 

MARSHMALLOW 

3.99 

4.0 

4.5 

CLEARWATER 

4.71 

4.9 

5.3 

HANDCAR 

4.19 

4.7 

(Dolomite) 

HARDHAT 

4. I5(mt8) 

4.6 

4.9 

(Granite) 

SHOAL 

4.62 

4.8 

4.9 

(Granite) 

DANNY  BOY 

3.06 

m 

3.6 

(Basalt) 

GNOME 

4.27 

- 

4.6 

(Salt) 

CLIMAX 

3.73 

• 

4.4 

(Granite,  Chemical) 

CHASE  III 

4.7  («>£■) 

• 

• 

(Water) 
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RELATIVE  MAGNITUC  P.  USE  OF  DATA  OF  FIXED  STATIONS 


Another  approach  to  establishing  relative  magnitudes  of  NTS  events  is 
to  compare  A/T  values  of  the  events  as  recorded  at  fixed  stations. 
Relative  magnitude  should  be  determinable  to  within  several  percent  by 
this  approach  and  one  measure  of  the  meaningfulness  of  the  overall 
scheme  here  employed  would  be  the  agreement  of  relative  magnitudes  of 
NTS  events  as  computed  by  the  two  methods  of  analysis.  Several  events 
are  analyzed  by  relative  amplitudes  at  fixed  stations  in  Table  IX,  all 
stations  being  restricted  to  ones  receiving  P^g  in  WUS  from  NTS 
events. 


TABLE  IX 

Comparison  of  m^i  and  Individual  Station  Relative  Magnitudes 


Event 

mb. 

Relative 
(CHIN  II) 

i  Rel 

MNNV 

stive  m 
WINV 

(CHIN 

KNUT 

II)* 

FSAZ 

Av 

6  m 

CHINCHILLA  II  2.73 

2.73* 

2.62* 

2.35* 

3.12* 

3.01* 

2.78* 

*.05 

KAWEAH 

3.15 

.42 

.73 

1.18 

-.10 

.03 

.51 

*.09 

STOAT 

3.38 

.65 

.89 

•• 

.59 

.62 

.70 

*.05 

CODSAW 

3.40 

.67 

.84 

1.10 

.34 

.19 

.62 

-.05 

FISHER 

3.45 

.72 

1.20 

.80 

.42 

.59 

.75 

*.03 

ACUSHI 

3.51 

.78 

1.11 

.96 

.35 

.84 

.82 

*.04 

DORMOUSE 

3.64 

.91 

1.30 

• 

.53 

.70 

.83 

-.08 

CASSELMAN 

3.64 

.91 

1 . 16 

1.49 

.61 

.68 

.99 

*.08 

PASSAIC 

3.69 

.96 

1.23 

1.34 

.72 

.75 

1.01 

*.05 

MADISON 

3.71 

.98 

1.08 

1.32 

1.07 

1.09 

1.14 

♦.  16 

♦CHINCHILLA 

II  data 

are  basis 

of 

comparii 

son  for 

other 

events. 
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In  Table  IX,  the  data  of  CHINCHILLA  II  are  the  basis  of  comparison.  Thus, 
in  the  row  tor  CHINCHILLA  II,  values  are  magnitudes  as  estimated  by  017.9 
while  in  the  rows  fot  other  events,  values  are  differences  bet’  ecn  magni¬ 
tude  for  these  events  and  CHINCHILLA  II  at  the  indicated  station  or 
stations,  based  on  measured  A/T  values.  Note  that  such  an  approach  does 
not  necessarily  eliminate  a  fu\l  magnitude  spread  in  the  data  (KAWEAH, 
for  example).  Column  8  is  an  average  of  the  relative  magnitude  values 
obtained  at  the  four  fixed  stations  (KNUT,  MNNV,  FSAZ,  and  WINV)  while 
Column  9  is  the  difference  between  this  relative  magnitude  estimate  and 
that  of  the  magnitude  values  of  this  report  (m^i).  In  only  one  case, 
MADISON,  does  a  value  in  Column  9  exceed  0.1  magnitude  unit.  The  conclu¬ 
sion  is  that  the  procedural  scheme  used  of  averaging  values  computed  by 
various  phase  formulas  has  passed  another  significant  test. 

APPLICATION  OF  THE  PHASE  CURVES  TO  EARTHQUAKES  OF  WESTERN  UNITED  STATES 

1.  General .  Problems  of  two  sorts  expectedly  arise  in  the  application 
of  these  tormulas  to  data  from  earthquakes.  First,  the  factor  of  rela¬ 
tive  excitation  of  phases,  significant  in  such  a  small  region  as  NTS, 
expectedly  introduces  inconsistencies  in  predicted  magnitudes  of  earth¬ 
quakes.  Second,  the  certain  definition  of  phase  for  earh  arrival 
measured  becomes  more  difficult  than  tor  NTS  shots.  For  shots,  phase 
definition  is  comparatively  simple  since  we  know  accurately  the  time  of 
origin  and  epicenter  and  2-second  differentials  in  travel-times  are  clear 
and  meaningful.  However,  the  programs  used  in  locating  earthquakes  do 
not  consider  the  presence  of  two  dominant  refracted  phases  in  the  dis¬ 
tance  range  300-1800  kilometers.  Ihe  attempt  to  locate  epicenters  on  the 
basis  of  a  formula  that  obscures  these  two  arrivals  by  treating  them  as 
one  or  as  an  invalid  two  can  lead  only  to  phase  confusion  and  lack  of 
definition  on  a  T  -  A/8.1  plot.  SHOAL  and  the  Fallon  Earthquake  form  an 
excellent  demonstration  pair  of  this  phenomenon.  These  r wo  events 
occurred  very  near  each  other  and  many  of  the  same  stations  recorded 
them.  The  T  -  A/8.1  plot  for  SHOAL  (see  shot  report)  clearly  shows  the 
several  refracted  phases.  Time  of  origin  and  epicenter  were  very 
accurately  known  and  were  rot  derived  from  the  travel-time  data.  The 
Fallon  Earthquake,  however ,  of  necessity  was  located  by  use  of  the 
travel-time  data  and  by  use  of  an  earth  model  which  did  not  represent  the 
velocity  structure  of  the  region.  Thus,  appreciable  uncertainties  in 
the  time  of  origin  and  epicenter  resulted  in  a  confused  and  virtually 
undecipherable  T  -  A/8.1  plot.  In  such  a  situation,  one  must,  almost 
of  necessity,  establish  arbitrary  ranges  for  the  applicability  of  each 
phase  magnitude  formula.  We  adopt  the  convention  that  mj  ^  applies  to 
all  stations  with  <  A  £  1000  filometeis,  to  all  stations  with 

1000  <  A  *  1800,  and  mi 0.3  and  mts  beyond  1800  kilometers.  The  breaking 
of  this  rule  occurs  when  obvious  reason  to  do  so  exists  and  we  note  such 
exceptions.  Expectedly,  the  shadow-zone  effects  of  the  Sierra  Nevada 
influence  relative  amplitudes  of  these  phases.  With  these  limitations 
in  mind,  we  selected  for  study  a  set  of  WUS  earthquakes.  Several  events 
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reported  by  the  L'S  Coast  and  Geodetic  Survey  (USC&CS)  to  he  of  magnitudo  5 
or  greater  augment  earthquakes  previously  reported  by  VELA  in  shot  report 
style.  As  an  interesting  sidelight  of  this  investigation,  note  that 
PSCSC.S  or  Pasadena  often  greatly  ovi rest imat e  magnitudes  of  small  earth¬ 
quakes  (5  or*  less)  relative  to  teleseismic  magnitudes.  Much  of  i  in i  •. 
clearly  duo  to  use  of  invalid  formulas  rr:  lat  ing  tuagn  i*-udc  ,  distance, 
and  A/I. 


Earthquakes  Studied. 

The  earthquakes 

investigated  arc  those  ot : 

> 

Feb 

1962  at 

18° 

12' 

'  N  104° 

54'  W 

5 

Feb 

1962 

38 

12 

•  107 

36 

Colons  Earthquake  Report 

15 

Feb 

1962 

36 

51 

112 

26 

r> 

Ft  b 

1962 

37 

04 

11? 

58 

2.') 

Feb 

1962 

45 

12 

U  1 

12 

hebgen  Earthquake  Report 

5 

Apr 

1962 

38 

24 

1 19 

18 

Bridgeport  Earthquake  Report 

14 

Apr 

1962 

41 

125 

27 

May 

1962 

31 . 

7 

115. 

6 

20 

Jul 

1962 
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values.  The  general  phenomenon  of  predicting  too  high  a  magnitude  by  use 
of  msr  at  less  than  1700  kilometers  becomes  as  apparent  as  for  NTS  shots. 
Fallon,  in  particular,  interests  us  as  its  pattern  of  computed  magnitude 
values  is  so  similar  to  that  of  SHOAL,  implying  essentially  identical 
energy  partition  functions  of  these  two  events  as  regards  P-wave  energy. 
The  events  ot  l A  April  1962,  27  May  1962,  6  November  1962,  and  24  October 
1964  display  excellent  adjustment  of  the  data  by  use  of  the  phase 
formulas  in  contrast  to  the  poor  adjustment  when  using  msr.  Events  of 
2  februa^y  1962,  21  August  1962,  !6  September  1962,  22  March  1964, 

31  March  1 964  and  July  :96c  gi  c  lower  *.*vc  males  of  magnitude  at  less 
than  approximately  L400  kilometers  than  at  greater  epicentral  distances. 
Note  that  these  events  occurred  either  immediately  west  of  the  Sierra 
Nevada  or  in  Mexico,  or  near  Queen  Charlotte  Islands. 
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3.  Small  Earthquakes.  The  quakes  °f  5,  15,  and  2S  February  1962, 

5  April  1962,  and  6  August  1963,  i.  e.„  Colona,  Kanab  aftershock,  Hebgen, 
Bridgeport,  and  o~e  other,  were  not  recorded  at  teleseismic  distances 
because  of  their  rmall  sire.  The  first  four  were  used  by  Mansfield  in 
his  discussion  of  AR.  Note  that,  of  all  of  them,  the  use  of  msr-  values 
predicts  a  teleseismic  magnitude  completely  inconsistent  with  the  short 
distant*  «  ovc  t  which  '■hese  events  were  recorded.  Adjustment  by  the  phase 
formulas  brings  the  individual  magnitude  data  into  more  reasonable 
internal  ugt eenu-nt  and,  in  addition,  yields  reasonable  magnitudes  in 
light  o(  the  short  rar.ge  of  recorded  signals.  A  more  detailed  discussion 
from  a  regional  pattern  point  of  view  seems  unwarranted  until  a  large 
sample  oi  nras  Iva  been  analyzed.  It  seems  obvious,  however,  that 
treatment  o(  the  A/T  data  of  earthquakes  of  WUS  by  the  WUS  phase  formulas 
yields  a  more  internally  consistent  body  of  data  than  that  obtained  by 
use  ot  msr  values. 

A,  Libya-  Har  T.jiu^e  oi  l 1  Feb  -  v  ai  y  1963.  As  a  !  rdl  note  on  this  point 

the  two  Lioyan  earthquakes  (figure  99)  studied  by  Gordon^  suggest  the 

appropriateness  of  application  of  these  formulas  (or  at  least  mg.5)  to 

the  Mediterranean  area.  For  the  two  events  reported  in  the  above- 

c  o  1 

mentioned  paper,  ATU  and  1ST,  S.8  and  10.6  ,  respectively,  yielded 
excessive  magnitude  estimates  by  use  of  msr  values  but  yield  values 
consistent  with  mt  s  by  for  1ST  and  ra^g  lor  ATU. 

USE  OF  P  PHASES  AS  DIAGNOSTIC  AIDS 

The  agreement  o»  m^Q  (s*e  Appendix),  myq,  mg^,  mjQ  5,  and  mLs  values 
for  both  explosior  s  a'd  earthquakes  means  that  relative  excitation  of 
these  phases  .  anno*,  b*  used  as  a  diagnostic  criterion  for  distinguishing 
explosions  and  *-'«irt h \  akr>s.  Li»»-le  dorbt  oists  that  patterns  of 
variable  relative  e«:itation  such  as  observed  for  BILBY  and  HARD MAT  will 
occur  m  *.-ar  thquakes  also  tut  a  pattern  separating  earthquakes  and 
explosions  dot?  not  exist. 

VELOCITY  SIRUCIUHE  0^  CRUST  OR  CRUST  AND  UPPFR  MANTLE  IN  EASTERN  UNITED 
STAIFS  AND  WF STERN  UNITED  STATES 

1.  EUS 

a.  limited  Pg  data  indicate  that  in  EUS  this  phase  is  never  the 
clearly  marked  phase  that  it  is  in  WUS.  The  events  picked  as  Pg  appear 
to  have  a  zero-distance  time- intercept  of  several  seconds  and  a  velocity 
of  about  6.4  kilometers  per  second.  However,  the  data  are  exceedingly 
poor. 

b.  Figure  23  is  the  basis  of  the  following  discussion,  being  based 
upon  numerous  shot  reports.  The  Ps.5  event  records  as  first  arrival 

to  2100  kilometers.  The  SS  Village  profile  (recorded  by  land  stations. 


^Gordon,  David  W, ,  1963,  Libyan  Earthquake  of  21  February  1963, 
USC&GS  Publ  i  *at  ion. 
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the  neatest  being  a  tcw  hundred  kilometers  from  the  shot  point)  has  an 
intercept  of  +9  seconds.  This  represents  a  path  which  includes  one  trip 
through  the  EUS  crust,  one  trip  through  an  oceanic  crust,  and  a 
horizontal  leg  of  a  few  hundred  kilometers  to  continental  structures  at 
about  8.1  kilometers  per  second  rather  than  the  continental  velocity. 
CNOME-East  shows  an  intercept  of  -*-11  seconds,  this  intercept  being 
representative  ot  two  trips  through  the  EUS  crust.  We  will  assume  that 
5.5  seconds  characterize  crustal  thickness  in  EUS,  the  3.5  second 
difference  of  the  SS  Village  Intercept  being  explained  by  the  two  last 
path  segpients  mentioned  above. 

c.  Very  simple  crustal  models  suffice  for  the  discussion  which 
follows.  The  discussion  illustrates  the  contrast  between  EUS  and  WUS, 
a  contrast  so  gross  that  slight  details  of  crustal  model  are  of  second- 
order  significance.  Therefore,  we  assume  a  velocity  of  6.0  kilometers 
per  second  above  the  Pg.5  refractor  and  obtain  a  thickness  of  46  kilo¬ 
meters  for  the  crust  in  EUS.  This  thickness  varies  as  other 
investigations  have  shown,  but  an  average  figure  suffices  for  this 
discussion. 


2.  WUS.  A  cursory  glance  at  Figure  23  reveals  immediately  the  gross 
contrast  between  travel-time  curves  for  EUS  and  WUS  in  the  range  150- 
2100  kilometers.  As  we  have  already  illustrated,  marked  differences 
occur  within  WUS  itself. 

a.  Pg  is  well  observed  in  WUS  out  to  1000  kilometers  and  has  a 
velocity  of  very  nearly  6.0  kilometers  per  second. 

b.  P79  is  widely  observed  in  WUS  for  NTS  events.  Its  intercept 
does  not  depend  upon  azimuth  from  NTS  and  is  about  +5  seconds.  Use 
of  a  layer  velocity  of  6.0  kilometers  per  second  indicates  a  depth  of 
26  kilometers  to  this  refractor. 

c.  The  P8.5  ■  NTS/Dakotas  profile  has  an  intercept  of  +14  seconds, 
representing  the  sum  of  1-way  trips  through  WUS  and  EUS  crustal  types, 
all  times  having  been  measured  in  EUS  for  this  phase  on  this  profile. 

If  5^  seconds  represent  a  1-way  trip  through  the  EUS  crust,  8^  seconds 
characterize  a  1-way  trip  through  the  WUS  crust. 

d.  For  P8.5  -  NTS/Washington,  an  observed  intercept  of  about 
■♦■18  seconds  should  characterize  a  2-way  trip  through  the  WUS  crust. 

A  value  of  9  seconds  for  a  l-way  journey  agrees  nicely  with  the  value 
of  8^  found  on  the  NTS/Dakotas  profile. 

e.  If  we  assume  an  average  value  of  8.0  kilometers  per  second 
for  the  velocity  in  the  layer  between  the  P79  and  Pg  5  refractors 

in  WUS,  a  thickness  of  approximately  130  kilometers  is  found  for  this. 
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layer.  The  Pg  tj  refractoc  in  WUS  then  occurs  at  a  dcptn  of  about 
150  kilometers  while  in  EUS  it  occurs  at  only  approximately  50  kilo¬ 
meters.  The  quest  on  immediately  arises  as  to  what  to  call  the  crust 
in  these  areas.  No  logical  basis  exists  for  calling  the  P7  9  refractor 
base  of  crust  in  WUS  and  the  Pg#5  refractor  the  equivalent  level  in  EUS. 
This  customary  approach  may  obscure  a  fundamental  and  profoundly  signi¬ 
ficant  contrast  between  these  two  parts  of  the  United  States,  No  logical 
explanation  of  the  supposed  thinning  of  the  crust  from  east  to  west 
exists.  The  general  conclusion  of  compensation  below  the  crust  has  been 
invoked  in  order  to  explain  the  conflict  between  gravity  and  seismic 
data.  Real  understanding  of  the  relations  between  these  two  types  of 
data  may  result  from  the  realization  that  the  same  velocity  surface  is 
present  in  WUS  as  is  termed  base  of  crust  in  EUS  but  is  approximately 
100  kilometers  deeper  in  WUS,  Probably,  calling  the  Pg  refractor  in 
both  areas  base  of  crust  would  constitute  a  geophvi-  ical  lv  meaningful 
definition. 

f.  P0  1  -  NTS/Northern  New  Mexico.  On  large  events  at  NTS,  a 
clear  8.0  -  8.1  arrival  is  recorded  as  f .ir  as  RINM,  the  intercept 
being  approximately  7.5  seconds  Detailed  inspection  of  records 
obtained  along  this  profile  indicates  absence  of  a  7,7  79  refractor. 

Indefinite  recordings  of  Pg  5  have  been  measured  at  a  lew  hundred 
kilometers,  really  good  recordings  not  being  obtained  until  the  profile 
extends  into  EUS,  The  Pg  j  refractor  apparently  replaced  the  P7  9 
refractor  in  this  region  and  is  at  a  somewhat  greater  depth,  the  values 
given  above  yielding  a  depth  of  36  kilometer i,  The  Pg  5  refraction  leg 
observed  in  Texas  and  beyond  has  an  intercept  of  18  seconds..  In  light 
of  the  fact  that  the  GNOME  East  profile  traversed  the  sam*1  region  with 
an  intercept  of  only  11  seconds,  the  suggested  explanation  is  that  the 
Pg  5  energy  observed  along  this  profile  for  NTS  shots  is  the  energy 
that  was  following  the  Pg  ^  refractor  m  WUS 

APPLICATION  OF  IDEAS  OF  THIS  PAPER  TO  OTHER  AREAS 

1.  General .  The  functional  relation  betw^n  A/T  and  A  for  a  given 
region  quite  obviously  depends  upon  the  veloctiy  structure  of  the 
region,  A  simple  velocity  structure  such  as  that  of  EUS  will  result 
in  a  simple  dependence  of  A/T  upon  A.  Conversely,  the  existence  of 
multiple  refractors  will  result  in  a  complex  dependence  of  A/T  upon  A. 
The  general  suggestion  is  that  Shield  or  Cratonal  areas  of  low  mean 
elevation  will  display  a  simple  A/T  versus  A  relationship  while  regions 
of  high  mean  elevation  on  a  regional  basis  will  display  complex 
velocity  structures  and  a  complicated  dependence  of  A/ I  upon  A. 
Therefore,  an  investigation  of  the  velocity  structure  of  a  region 
will  lead  to  a  basic  understanding  of  the  operative  A/T  versus  A 
relationships.  Such  velocity  data  cannot  be  based  upon  an  epicentral 
location  and  origin  time  based  in  turn  upon  an  assumed  velocity 
structure.  Thus,  only  data  for  explosions  are  acieptable  for  such 
an  analysis. 
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2.  Europe.  The  Haslach  travel-time  curve  for  Europe  allows  an 
immediate  prediction  of  the  expected  A/T  versus  A  relationships  for 
this  region.  The  data  on  which  the  travel-time  curve  of  Figure  23 
is  based  were  obtained  as  the  result  of  an  explosion  near  Haslach 
in  Southwest  Germany  (48°16'N  08°07'E).  Note  that  a  single  phase 
with  uniform  apparent  velocity  was  recorded  from  200  kilometers  to 
2000  kilometers,  a  pattern  strikingly  similar  to  that  observed  in 
EUS.  It  seems  reasonable  that  A/T  will  be  found  proportional  to 
very  nearly  the  Inverse  second  power  of  A  in  Europe.  There  seems 
to  be  little  possibility  of  a  pattern  similar  to  that  of  WUS.  No 
shadow  zone  for  P  phases  will  be  found  in  this  area. 

3.  Asia.  The  USSR  has  not  released  adequate  explosion  data  to 
permit  a  prediction  of  the  expected  A/T  versus  A  relationships. 
Obtaining  such  data  is  important  as  it  would  greatly  advance  our 
understanding  of  the  problems  to  be  expected  from  observatories  at 
specific  localities  within  the  area  of  interest., 

CONCLUSIONS 

1.  Contrary  to  the  curve  presented  by  Richter**,  Pn  amplitudes  in  WUS 
In  the  range  1000-1700  kilometers  are  not  obeying  a  shadow  zone  inter¬ 
pretation  as  envisaged  by  Gutenberg.  There  is  no  evidence  for  a  broad 
interval  of  several  hundred  kilometers  in  which  P  amplitudes  are  very 
small.  There  is  a  large  jump  in  amplitude  at  approximately  1000  kilo¬ 
meters  due  to  Pg'^  becoming  the  first  arrival.  For  some  uncertain 
reason,  the  Pg  ^’refractor  and  associated  refracted  phase  were  never 
recognized  in  WUS  by  Gutenberg  and  the  high  amplitudes  of  this  phase 
were  either  missed  or  misinterpreted.  Signal  levels  decrease 
continuously  from  approximately  1000  kilometers  outwards  except  for 
one  or  more  comparatively  small  jumps  in  amplitude  associated  with 
the  beginning  of  the  first  arrival  of  deeper  refractions.  These  small 
steps  in  amplitude  (such  as  from  Pg  5  to  P^q  cannot  be  termed 
shadow  zone  phenomena  in  the  Gutenberg  sense  or  in  any  real  sense. 
Measurement  of  P  amplitudes  in  the  range  1000-2000  kilometers  are 
more  desirable  than  in  any  more  distant. 

2.  In  EUS,  amplitudes  increase  continuously  from  1000  kilometers  to 
the  epicenter.  There  is  total  absence  of  any  phenomenon  which  could 
be  interpreted  as  a  shadow  zone  for  P> 

3.  In  WUS,  the  amplitude  discontinuity  in  the  neighborhood  of  1000 
kilometers  is  abrupt  and  associated  with  the  first  arrival  of  a 
deeper  refraction  (P8,5)>  replacing  the  weaker  P7.9  event.  The 
P7.9  event  can  occasionally  be  followed  as  second  arrival  beyond 

A 

Richter,  loc.  clt.  and  Figure  100,  this  report. 
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the  point  where  Pg^  becomes  first  arrival.  The  P8.5  arrival  in 
WUS  appears  to  obey  essentially  the  same  decay  law  with  distance 
as  does  the  P8.5  refraction  in  EUS. 

4.  The  A/T  versus  A  table  used  in  the  VELA  shot  reports  is  grossly 
inaccurate  at  distances  out  to  15°-20°  for  virtually  all  events,  even 
for  WUS.  Errors  in  magnitude  estimates  at  single  stations  of  as  high 
as  l.S  magnitude  units  occur  when  this  lav  is  employed.  Extensive  con¬ 
firmation  of  this  failure  of  the  shot  report  table  or  for  any  single 
curve  for  both  explosions  and  earthquakes  is  contained  in  this  report. 

All  previously  studied  explosions  and  earthquakes  plus  several  additional 
earthquakes  have  been  recomputed  according  to  the  scheme  proposed  in 
this  paper.  Romney's  procedure  for  computing  m  for  small  events  by  use 
of  P7.9  data  was  calibrated  against  teleseismic  data  of  large  events  and 
thus  avoided  the  inaccuracies  of  the  shot  report  computations. 

5.  A/T  data  for  the  Pg  phase  in  WUS  are  as  useful  and  as  meaningful 
for  magnitude  computations  for  both  earthquakes  and  explosions  as  are 
the  refracted  P  phases  (see  Appendix). 

6.  Magnitude  computations  based  on  refracted  phases  can  be  made 
essentially  independent  of  distance  by  proper  calibration  of  the 
crust  and  by  establishing  the  correct  energy  partition  function 
for  the  several  phases  generated.  Such  energy  partition  functions 
appear  quite  variable  and  must  be  influenced  by  focal  and  near-focus 
conditions. 

7.  Epicenter  location  of  unknown  events  by  use  of  regional  and  near- 
regional  data  will  be  markedly  improved  only  when  crustal  models 
expressing  the  multiple  refractors  present  in  a  region  are  used  as 
the  basis  of  analysis.  Any  effort  to  attach  regional  corrections 
must  fail  in  general  application. 

8.  Analysis  of  the  amplitude/period  data  of  any  region  in  the  range 
200-2100  kilometers  is  only  meaningful  after  elucidation  of  crustal 
structure  and  establishment  of  the  general  energy  partition  function 
and  decay  rates  of  each  phase.  It  appears  probable  that  velocity 
structure  will  characterize  decay  rates.  Thus,  the  Haslach  travel¬ 
time  curve  for  Europe,  based  upon  surface  explosion  data,  is  very 
similar  to  that  of  EUS.  The  implication  is  that  mg, 5  may  be  appli¬ 
cable  to  A/T  data  of  European  seismic  events  in  the  regional  and 
near-regional  distances  but  only  investigation  of  A/T  data  of  Europe 
will  indicate  the  actual  relationships.  The  situation  in  the  USSR 
is  not  certain  and  will  not  be  clarified  until  data  obtained  by 
observing  Russian  explosions  become  available. 
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APPENDIX:  ANALYSIS  OF  Pg  DATA 

1.  General .  P6.0  (Pg)  was  observed  widely  and  well  throughout  WUS  to 

distances  of  1000  kilometers  with  a  velocity  of  6.0  kilometers  per 
second.  Pg  gives  rise  to  higher  A/T  values  than  other  P  phafses  in  WUS. 
In  many  cases  it  can  be  recognized  by  waveform  character  even  though, 
at  other  than  very  near  stations,  it  is  a  second  arrival.  In  WUS  it 
decays  as  the  third  power  of  the  epicentral  distance  and  its  velocity 
is  very  near  6.0  kilometers  per  second.  It  does  not  occur  beyond  WUS 
and  loses  definition  beyond  about  1200  kilometers  in  any  azimuth  from 
NTS.  One  wishes  to  know  whether  one  can  interpret  Pg  A/T  data  in  a 
magnitude  sense  to  the  extent  that  the  Pn  phases  have  been  successfully 
interpreted.  If  so,  two  useful  and  significant  results  accrue. 

Firstly,  one  can  obtain  confirmation  of  our  preposed  A/T  versus  A 
versus  m^i  interpretive  scheme  from  an  independent  phase.  Secondly, 
such  a  uniformity  of  behavior  would  imply  one  ;ould  use  Pg  for  estima¬ 
tion  of  magnitude  of  earthquakes  with  as  much  or  more  reliability  than 
the  Pn  phases  as  the  scheme  of  identification  of  Pn  phases  for  an 
unknown  epicenter  may  prove  difficult  while  one  can  generally  identify 
Pg  by  record  character.  If  the  beginning  of  P  is  measurable  to  within 
two  or  three  seconds,  better  epicenters  probably  could  be  obtained  by 
use  of  this  phase  than  by  use  of  the  Pn  phases  in  WUS.  Pn  phases  as 
first  arrivals  are  timed  to  a  fraction  of  a  second  but  the  proper 
classification  of  the  phase  measured  may  be  so  uncertain  and  the 
crustal  velocity  structure  assumed  so  inaccurate  that  one  can  introduce 
errors  of  several  seconds  into  the  calculations  of  epicenters.  When 
using  Pg  in  WUS  an  extremely  simple  but  valid  model  can  be  used  for 
computation  purposes. 

2.  Analysis  of  Pe  Data.  The  procedure  followed  here  begins  by 
estinating  the  best  value  of  A/T  at  500  kilometers  from  the  A/T 
versus  A  plots  for  Pg  in  the  shot  reports.  See  Figures  103  through 
106  for  examples  of  A/T  versus  A  plots  of  Pg  data.  We  plot  these 
(A/T)500  values  against  radiochemical  yield  in  Figure  107  and  against 
m^t  values  of  this  report  in  Figure  108.  Compare  Figures  107  and  108 
to  see  that  the  (A/T)500  values  of  Pg  display  a  linear  relationship  to 
the  magnitude  values  of  this  reports  being  independent  of  shot  medium. 
These  two  figures  confirm  the  usefulness  of  Pg  amplitude  data  for 
estimating  magnitude  and  support  the  relative  magnitude  determinations 
of  this  report.  The  resulting  P^q  equations  are: 

m6.0  "  1.28  ♦  1.14  log  (A/T)500  (13) 

m6.0  “  “8-00  *  lo8  A/T  ♦  3.42  log  A  (14) 
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The  curves  of  Figure  100  for  P7  tg  can  be  used  for  Pg^  if  magnitude 
values  are  decreased  by  1.0  units.  In  other  words »  P,  ,  gives  A/T 
values  10  times  those  of  P79  at  the  same  distance. 

3.  Pg  Data  for  WUS  Earthquakes.  A/T  data  for*  the  P„  phase  are 
available  for  several  earthquakes  of  WUS.  The  computed  mgQ  values 
follow: 


Earthquakes 

(A/T) 500 

V 

m6.0 

Colons  (5  Feb  1962) 

40 

3.1 

3.11 

Hebgen  (25  Feb  1962) 

50 

3.2 

3.22 

Cache  Creek  (30  Aug  1962) 

2500* 

4.3 

5.16 

Cache  Creek  A/S  (5  Sep  1962) 

500* 

4.1 

4.36 

Bridgeport  (5  Apr  1962) 

42 

3.03 

3.29 

Fallon  (20  Jul  1962) 

240 

4.4 

4.16 

*Based  on  WUS  data. 


The  Pg  data  demonstrate  their  usefulness  for  the  estimation  of  mag¬ 
nitude  of  WUS  earthquakes.  Remember  that  we  derived  this  scheme  for 
Pg  by  calibration  against  mt  values  of  P  phases  and  noting  the 
linear  relationship  between  log  Y  and  log  (A/TJjqq.  Therefore,  the 
agreement  at  magnitudes  above  4  occurs  independently  of  Pn  values 
inside  2000  kilometers  and  the  agreement  at  lower  magnitudes  depends 
upon  control  other  than  previously  computed  magnitude  values. 
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